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o 

' For understanding the origin of anisotropies in the cosmic microwave background, rules to 

' construct a quantized universe is proposed based on the dynamical triangulation method of 

' the simplicial quantum gravity. A d-dimensional universe having the topology D'^ is created 

' numerically in terms of a simplicial manifold with d-simplices as the building blocks. The 

' space coordinates of a universe are identified on the boundary surface S'^~^, and the time 
coordinate is defined along the direction perpendicular to S'^^^. Numerical simulations are 

^ ' made mainly for 2-dimensional universes, and analyzed to examine appropriateness of the 

' construction rules by comparing to analytic results of the matrix model and the Liouville 



o 



o 



theory. Furthermore, a simulation in 4-dimension is made, and the result suggests an ability 
to analyze the observations on anisotropies by comparing to the scalar curvature correlation 
of a S'^-surface formed as the last scattering surface in the universe. 

KEYWORDS: Origin of Universe, CIVIB, Inflation, Quantum Gravity, Dynamical Triangulation, 
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■ 1. Introduction 

5h ' Recent observations of anisotropies in the cosmic microwave background (CMB) by COBE 

and WMAP^"^ extend our knowledge on the history of the universe beyond the big bang. They 
imply as follows: about 14 billion years ago, the universe started expanding with extremely 
rapid acceleration, now known as inflation,'^'^ from a micro-universe of quantum mechanical 
scale. Because of the expansion speed exceeding the light velocity, the pattern of original 
quantum fluctuations in space could not be disturbed dynamically, and was kept intact dur- 
ing the expansion something like a picture printed on an inflating balloon. After a little 
while(10~^'^ ~ 10^^^ sec) the expansion was slowed down through coupling to matter fields, 
and fluctuations developed in the space curvature transferred to matter density fluctuations. 
They eventually are observed as the temperature anisotropies in the CMB, caused by the 
photon energy dissipation in fluctuating media known as the Sachs- Wolfe mechanism.^ 

Although the observed data are rather well-understood by known physical processes with 
several phenomenological parameters, questions such as how and what kind of fluctuations 
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were created initially, and why the universe could expand so rapidly arc left to be the matters 
which essentially belong to the fundamental lows of space-time. In the history of studies on 
quantum gravity it is the most remarkable occasion where the theory is required to be phe- 
nomenological, apart from describing quantum fluctuation of geometry, the graviton, whose 
appearance is beyond the current experimental scope. It is a common belief that there exists 
no consistent theory that combines quantum mechanics and gravity at the moment. Among 
various efforts to establish quantum theory of gravity the super-string is regarded to be the 
most promising candidate as far as quantization of the graviton is concerned. It is extraordi- 
nary that the superstring theory can unify all the known interactions on the way of quantizing 
graviton. However, it does not say much about how the space-time is created. While absence 
of a consistent theory of quantum gravity is true in 4-dimension, we certainly have consistent 
theories of quantum gravity at least in 2-dimension. One of them is a standard field theory 
called the Liouville field theory,^ and the other is the matrix model. ^ In addition, the latter 
model has a numerical realization known as the dynamical triangulation (DT).^ 

The existence of analytic theories is closely tied to the dimensionality of space, and it is 
not a straightforward task at all to extend those theories into higher dimensions. The con- 
formal gravity-'^'^'^^ is one of possible generalizations of the Liouville theory into 4-dimension. 
It is free from space singularities such as the black hole of quantum mechanical size, and 
claimed to be renormalizable. However, it is unfortunate that the non-unitary nature inher- 
ited in any higher derivative field theory masks charms of the theory. ""^^ An extension of the 
discretizcd space model into higher dimensions beyond the matrix model has been proposed 
for the 3-dinicnsional space by adopting the 6j-symbols as the elementary unit of space. ""^^ 
In an analogous manner it may be possible to extend further to the 4-dimensional space by 
introducing an appropriate elementary unit.^^ The simplicial quantum gravity is one of the 
constructive generalizations in the discretized space approach, selecting a d-simplex as the 
elementary block of the d-dimensional space. Numerical simulations have been performed for 
the 2, 3, and 4-dimensional simplicial spaces, and the grand canonical Monte Carlo simu- 
lations have produced significant results not only for the 2-dimensional case,^^ but also for 
4-dimension^^ in calculations of the string susceptibility which plays the role of a universal 
observable of the space geometry. As far as the possibility of extending the method into 4- 
dimension is concerned the simplicial quantum gravity is one of the promising cases so far 
exhibiting its applicability. 

In modern field theories, quantization is often carried out by the Feynman path integral 
in which the Euclidean metric is adopted for definiteness of integrations. For a system having 
the stable vacuum it is always possible to transform the metric from the Euclidean to the 
Lorentzian and vice versa, by picking up one coordinate, which is regarded as the time coordi- 
nate, and making the Wick rotation. However, this procedure cannot be applied for a system 
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without the stable vaeuum. Quantum gravity is one of the systems having non-stationary 
vacuum of expanding space. Even though there exists no standard method of quantizing such 
a system, we employ the path integral method as a possibility to a space with open topol- 
ogy. While our final goal is to establish a method to construct the quantized universe in 
the 4-dimensional space-time, the following discussions are mainly proceeded on two dimen- 
sional space concerning the difficult circumstance of knowing no guiding theory. Besides its 
simplicity, analytic studies compiled extensively in 2-dimension arc inevitable for achieving 
some hints required in this totally new subject. We would like to stress, however, that the 
constructive nature of the method will make it readily applicable to higher dimensional cases 
as well. 

In the next section, we review the DT method mostly concentrating on the 2-dimensional 
space with the 5^ topology, and show that it provides a natural extension to a system with the 
open topology D"^. In Sec. 3, we shall practice to build an ensemble of 2-dimensional quantum 
universes based on two fundamental rules. In order to check the construction method we 
calculate observables of the space for systems with the stable vacua, by tuning the cosmological 
constant so that we can compare them to analytical results of the matrix model and the 
boundary Liouville theory. -"^^ In Sec. 4 we define the time and space coordinates of fluctuating 
manifold, and calculate the two point correlation function for the expanding universe. In 
the last section we give a summary and discussions on future directions such as a method 
extending into 4-dimcnsion and problems related to the introduction of matter contents of 
the Tiniverse. There, we show a numerical result of the sinuilation in 4-dinicnsion, which clearly 
indicates the applicability of the present method to the CMB anisotropics. 

2. Methods of simplicial quantum gravity and their extension 

The dynamical triangulation (DT) has been a standard method of non-perturbative quan- 
tum gravity^^ known as the simplicial quantum gravity. It was intensively studied in 2- 
dimension mainly because of its close connection to the string theory. Although we cannot 
expect new results from the 2-dimensional DT more than those already known from the ma- 
trix model and the Liouville theory, the achieved consistent results encourage us to adopt this 
method as a reliable technique for investigations beyond 2-dimension.^^ One of the drawbacks 
of the simplicial quantum gravity used to be the restriction to spaces with closed topolo- 
gies with the Euclidean metric. In the foUowings, we would like to show how to extend the 
d-dimensional DT technique to a space with the open topology D'^ for d > 2. Besides the 
space-time topology the dynamics requires a proper signature of the metric, which is known 
to be Lorentzian in this universe. We regard the selection of the signature should come out nat- 
urally when we consider dynamics, since the correct signature, which is required by dynamical 
laws, supposed to emerge in the process of establishing the space-time. 
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Quantum gravity is formally defined in a form of the path integral, 

Z = J V[g,,]eM-S[9ab])- (1) 

Here, we consider the simplest case that the metric tensor gab is the unique dynamical variable. 
The space is assumed to have the dimension d and the topology 5^ unless specified. The 
simplicial quantum gravity interprets the path integral to be a sum over all the distinct 
configurations of space discretized by d-simplices. For a d-dimensional manifold the action 
S[gab] is given by a linear combination of the variables, {Nn} with n = ■ • • d. It is considered 
to be natural, since we regard the diffeomorphysm invariance as its independence of the local 
geometry, or in other words, local variables. Thus only global variables such as the total 
number of n-simplices Nn are proper, and linearity comes from the extensive property of the 
action function. The (d + 1) variables {iV„} are not all independent. The Euler relationship 
and the manifold conditions reduce the number of independent variables. For example in 2- 
dimension, they are No — Ni + N2 = x with x being the Euler characteristic, and 3N2 = 
2Ni. This makes the action to be S^"^^ = IJLN2 with a parameter |Lt, which we call the lattice 
cosmological constant. 

The set of distinct triangulations is collected computationally by the Metropolis Monte 
Carlo method creating a Markov chain of configurations {oj}, 

ao ^ • • • ^ ^ Oj+i • • • . (2) 

Among possible methods of generating configurations one after another, we employ the so- 
called (p, g)-moves which allow us a systematic extension into higher dimensions. These moves 
are known to conserve the topology. There arc three types of (p, g)-moves in 2-dimension as 
shown in Fig. 2. From the figures it is clear that p and q stand for the numbers of triangles dis- 
appeared and appeared by the (p, g)-move, respectively. The Metropolis Monte Carlo method 
employs the importance sampling, accepting or non-accepting a new configuration created 
from an old configuration by a randomly selected move according to the detailed balance 
condition, 

Wa^lS = W/s^a, (3) 

Ha Up 

with the given probability pa for a configuration a. The factor n^, counts the number of 
possible moves starting from a configuration a. This is an important factor for simulations 
with varying the number of the degree of freedom in order to obtain the proper statistical 
weight for each configuration. 

The number of distinct triangulations is known to grow exponentially in N2 for 2- 
dimension, and the partition function behaves as 

Z~iVr'exp{-(/x-/ie)iV2}, (4) 
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Fig. 1. (p, g)-moves in the 2-dimensional simplicial gravity. 



where 7 is historically called the string susceptibility. As [i approaches to a critical value [ic 
from above, the average value of N2 tends to diverge as (// — //c)~^) and the system reaches 
to the continuum limit where the simulation meets the Liouville field theory with the field 
theoretical cosmological constant, 



parameter /i below /ic, configurations with more triangles become favorable, and rapidly 
increases as a simulation proceeds. Although it is beyond cases where the standard quanti- 
zation procedure is applicable, wc regard the partition function obtained by this method as 
the quantization of a system with an unstable vacuum, simply because the method is a trivial 
extension of quantization for a system with the stable vacuum. 

Another important setup we need to modify for a practical application of the method 
to a system with an unstable vacuum is the topological structure. The simplicial quantum 
gravity usually consider the space manifold within a closed topology. Of course we may regard 
it corresponds to a case of oscillating universe, if we take the space topology seriously as a 
real universe. However, it is widely believed that our universe has an open topology and it 
is expanding indefinitely in both space and time. Then, we wish to extend the DT method 
so that we can create the space with a boundary such as the d-disk D^. In order to show 
how it is possible, let us take the 2-dimensional simplicial gravity as an example. Among the 
three (p, g)-moves we examine the (l,3)-move, firstly. This move can be described either by 
dividing a triangle into three triangles by adding a vertex in the triangle, or by attaching 



A = 



(5) 




for space with a fixed volume A = ^a?N2, where a is the lattice constant. If we choose the 
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Fig. 2. {p,q)-moves in and [AF, A^J-moves in D^. 

one tetrahedron on a triangulated surface with one triangle face-to-face, which we call the 
[l,2]-moves representing changes of the volume AV and the surface area AS', as shown in 
Fig. 2 (a). 

In the similar manner the (2, 2)-move as well as the (3, l)-move can be done by attaching 
a tetrahedron with 2 or 3 triangles glued to those neighboring triangles of the S^-surface, 
which can be written in terms of [AV^, ASJ-moves as [1,0] and [1,-2] shown in Fig. 2 (b) 
and (c), respectively. In addition to these three moves it is also possible to make the three 
types of {p, g')-moves by detaching one tetrahedron from a S^-surface in just the opposite 
manners to attaching a tetrahedron.... It is clear these six moves constitute a set of ergodic 
moves of triangulations for the 3-dimensional space with topology. Each move changes the 
number of total tetrahedron, N3, by one unit, and at the same time it makes one of the three 
{p, g)-movcs on the boundary S^-surface. This method can be extended straightforwardly to 
other dimensions, such as the 2-dimensional disk, D^, and the 4-dimensional disk, as we 
will see in the following sections. Each case constitutes the four and eight [AV, ASJ-moves, 
respectively. 
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Fig. 3. 2-diniensional [AV, A5]-moves. The vertices and links on the boundary are denoted as 
white circles and the dash lines. 



3. Creation of simplicial universe 

When the lattice cosmological constant is below the critical value, the simplicial space 
tends to expand exponentially as a simulation proceeds. If we can regard this phenomenon 
to be the simplicial space analog of what happened in our universe at its very beginning, 
the inflation, the simulation can serve as a good numerical tool to study early stages of 
the evolution of universe. However, when we try to simulate the birth of a universe we find 
ourselves in a state of complete ignorance: we know neither dynamics nor initial condition. 
We have to recognize that there is no dynamical law to start with, since any physical law 
is considered to appear only at or after the space-time is settled in its classical sense. Also, 
there is no universe but a fluctuation at the beginning, and thus no initial condition. Under 
these circumstances we are forced to accept a priori elementary rules, but they should be 
as few as possible with the least prejudice. We shall employ an analogy of the creation of 
the natural number space where two rules known as the Peano axioms are essential: 1) there 
is an elementary unit, which is denoted by 1, and 2) there is the next neighbor S{a) of a 
natural number a. We are surprised by many profound mathematical laws such as Fermat's 
last theorem hidden in this simple space. 

The axioms for the d-dimensional simplicial universe we propose are constituted by 1) there 
is an elementary unit of universe, which we choose a d-simplex, and 2) there is a set of neighbor 
universes -S'(q;), constructed by either attaching or removing an elementary d-simplex to or 
from a universe a. There will be many neighbor universes possible. Among them we accept 
only those universes which fulfill the manifold conditions: 1) at most two d-simplices can 
attach through one {d — l)-simplex called the face, and 2) d-simplices sharing one vertex form 
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Fig. 4. 2-dimensional universes up to 4 triangles 

a (d+ l)-disk, or a {d+ l)-semi disk for the case when the vertex is on the boundary. According 
to our least prejudice principle this kind of manifold conditions may not be necessary, but we 
expect those universes which can grow indefinitely large will naturally select configurations 
following appropriate manifold conditions from mathematical consistency. We shall impose 
these manifold conditions at the beginning just for computational simplicity. However, the 
condition may not be unique. In Appendix B we define a different type of moves which create 
a manifold having more severely restricted configurations, while the manifold belongs in the 
same universality class. There even exist manifold conditions which allow topology changes. 
Inclusion of various possible topologies in the partition function is beyond the scope of this 
paper, and we restrict moves which conserve the topology. 

Now, let us start making a two dimensional universe as an exercise. The 2-dimensional 
minimum universe is a triangle, which acts also as the elementary unit. There exist four 
moves satisfying the manifold conditions as shown in Fig. 3, where we have denoted each 
move as [AN2, ANi] with the boundary link number Ni, instead of Ni, which is related as 
Ni = 2Ni — 3N2. Quantization of the simplicial universe is carried out by constructing the 
partition function to be a set of all the possible distinct triangulations. Possible configurations 
are collected as a Markov chain formed by selecting one of the four moves one after another, 
with appropriate weights imposed by the detailed balance condition. In Fig. 4, we show first 
few possible configurations up to universes made of 4-triangles. A configuration is specified 
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Fig. 5. Number of configurations for 
fixed I — 3. 



Fig. 6. Number of configurations for 
fixed A = 20. 



by four quantum numbers, {N2, Ni, k, s). The third quantum number k represents the type 
of a diagram. For example, they are assigned to be (a) and (b) in Fig. 4 for universes with 
= 4, A'^i = 4). The forth letter s specifies the point group quantum number of a diagram 
with multiplicity m^, which is A'^i for most of the cases except point symmetric diagrams such 
as Fig. 4 (b) which has = 1. 

In general, moves are independent of the quantum number s, and when we denote a set of 
the first three quantum numbers as a = {N2, Ni, k}, the detailed balance condition reduces 
to 

Parria Pbrrib , . 

Wa^b = Wb^a, (6) 

Ua rib 

where is the point group symmetry factor of the configuration a. For configurations with 
large there are seldom diagrams having point symmetries, and we can approximate = 
A^i, which we assume in our simulation. This assumption is checked by comparing to the 
matrix model prediction (Appendix A) obtained from expansion coefficients of the generating 
function restricting diagrams without tadpoles and self-energy insertions, for a fixed peripheral 
length I in Fig. 5 and a fixed area A in Fig. 6. 

The distribution is obtained by counting the number of states with equal weight pa = 1 
for each distinct configuration. The equal probability for each distinct configuration is also 
a natural rule for the quantized universe according to our least prejudice principle. Besides, 
we may give an arbitral probability {pa} through an appropriate action function Sa as Pa = 
exp(— 5a) depending on the purpose of a simulation. (Here, it may not be relevant to call Sa 
the action for such a system without dynamics. Nevertheless, we adopt the name only because 
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Tri-critical Point ^^=1. 1246(1), ^^^=0.837(1) 




Fig. 7. Three types of Universes. 

of the familiarity accustomed in statistical physics.) In d-dimension it is written as {Ni}, as 

d 

~^/i«iVi. (7) 

i=0 

By the topological and the manifold constraints, parameters in the action function will be 
reduced. For the disk topology only the Euler relation is required to satisfy and the action is 
expressed as 

5(2) =^iV2 + ^^iVi, (8) 

in 2-dimension, where we use A'^i instead of the total number of links A'^i. We remind these 
parameters are introduced only for the convenience to search possible types of universes in 
the followings. 

In the two dimensional parameter space {fi, the numerical simulation reveals two 
critical lines separating into three regions corresponding to three types of universes as shown 
in Fig. 7. When /x is greater than the critical values fic, the volume tends to decrease resulting 
a collapsing universe having a fate to disappear, while when fi and fi^ are smaller than the 
critical values both the volume and the surface area grow and the universe keeps expanding. 
The region where < /Xc and fj,^ > /x^ the volume increases while the surface decreases 
and the universe looks like an elongating tube with fluctuating radii. The critical parameters 
/Xc = 1.1247(1) and /if = 0.837(1) at the tri-critical point coincide to the values predicted 
by the matrix model, which we briefly describe in the Appendix A and B for two types of 
manifolds. 

In the numerical simulation, /Xc can be easily identified by inspecting the tendency of 
volume increase. However, the critical point /xf is not easily determined, because it does not 
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Fig. 8. The surface area distribution at the tri-critical point for a fixed volume A — 100, 1000, and 
2000. 



cause any divergence like /i. We fix this parameter so that the surface area distribution for a 
fixed volume coincides to the volume-area distribution function predicted by the matrix model 
as derived in Appendix A, 



3^ ) ~ ^'^^^^l exp{-(^^ 



Mf)iVi}. (9) 



When we identify the volume and the peripheral length of the boundary Liouville theory and 
the matrix model as ^ = ^a?N2 and / = aiVi, the distribution function (9) is rewritten as 



4 

2.5)0.5 



exp 



1 f 



exp(— \A) exp(-A^O 



(10) 



with = {ji^ — iJL^)/a, which differs from the Liouville theory prediction (C • 10) in two 
points, namely the sign in front of the exponent of /, and the factor 2 in the coefficient of ^. 
These differences is considered to be due to the way of counting the distinct configurations. In 
the boundary Liouville theory those configurations which differ only in numbering of external 
lines are thought to be equivalent,^*' and the exponent takes minus sign (Appendix C). While 
in the matrix model the point group symmetry factor is taken into account, and it chooses the 
plus sign (Appendix A, and B). The coefficient of ^ depends on the way to choose acceptable 
diagrams, and in facts, the Tuute algorithm gives 192/625 instead of 1/3 as described in 
Appendix A, and B. We can perform simulations either by the matrix model convention or 
the Liouville theory convention selecting the multiplicity factor in Eq. (6) to be A^i or 
1, respectively. In any way, this ambiguity will not cause any problem practically, since the 
multiplicity of a universe cannot be observed. 

The parameter h appearing in the distribution function is related to the background charge 
defined hy Q = h + in the Liouville theory. The exponent of A is given by — + |^ 
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and for the disk topology the Eular characteristic X = 1- For the pure gravity with the central 
charge cl = 26 it is known to be 6^ = 2/3. The numerical simulation shows the value very close 
to the theoretical prediction —2.5. The coefficient of in the exponent near the tri-critical 
point of the distribution function appears in the simulation as 0.30(4) as shown in Fig.'^ which 
is consistent to the theoretical value 1/3. From these agreements between simulations and 
theoretical predictions in the distribution function, we feel sure that the algorithm can indeed 
reproduce the universe predicted by the 2-dimensional quantum gravity. 

4. Introducing time and space 

Without knowing what is time the definition is inevitably heuristic expecting it will appear 
naturally reflecting the specific property of a universe. The reason why there are two types of 
directions in coordinates, time and space, is a reflection of the specific topology of our universe. 
For a universe with the disk topology D'^ it is natural to identify the space coordinate on the 
boundary S'^~^, and the time coordinate along the direction perpendicular to the surface. A 
definition of the physical time can be found in the trivial relationship between the volume 
and the area of a universe, 

Vit)= f S{t')dt', (11) 

where V{t) and S{t) represent the volume (~ N^a'^) and the surface area (~ N^-ia'^'^), 
respectively. In the numerical simulation for the expanding phase both the volume and the 
area increase as the simulation proceeds, which we record as functions of the diffusion time r, 
defined by counts of the number of throwing dice in a Monte Carlo simulation. Above relation 
can be rewritten in terms of r as 

S{r') dr' ^^^^ 
indicating the growing rate of physical time is proportional to the volume increasing rate per 
unit area. The physical time we define has the unique direction of arrow only in the expanding 
universe which occurs when the cosmological constant is below the critical value. It becomes 
negative for the case if the volume decreases in r. 

Simulations of the 2-dimensional universe in the expanding phase show the surface area 
increases approximately proportional to the volume as Fig. 9, 

V{t) ~ a-^S{T). (13) 

Thus, from the definition of the physical time, we obtain 

t = a-Hog{S{T)ISm, (14) 

or S{t) = S'(0)e'^*, exhibiting the exponential expansion of the universe. As a matter of fact 
the inflational expansion is a phenomenon we expect in the classical Liouville theory with the 
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Fig. 9. The ratio of V/ S for the expanding phase. 



negative cosmological constant as shown in the Appendix D. A homogeneous classical solution 
for the expanding mode is given by 

^2bm = a\^)e^'^'i>^^\ (15) 

where a(r/) = (cosw??)"^ is the conformal scale factor at the conformal time rj. The surface 
area of a homogeneous universe is the cross-section at r] given by L^/ cos(u;r/), where is the 
initial size along the conformal space coordinate ^. The volume of the universe at r/ is then 
given by 

-L^dT?'. (16) 



/ 



Iq cos^ {torj' 

We can obtain the physical time through our heuristic definition, using the volume and the 
surface area, which gives 



t 



1 



X 1 + sinw?? 

-log . (17) 

U) cos UJT] 

It coincides with the standard relation between the conformal time and the physical time, 
dt = a{rj)dr]. Once the time coordinate is determined the space coordinate is defined on 
the equal time hyper-surface, S"^~^. In the 2-dimensional universe it is unique to choose the 
coordinate axis along the periphery. 

The simplest definition of distance is the geodesic distance counting the number of links 
on the hyper-surface connecting two vertices with the shortest path. In terms of the distance 
we can measure the correlation function between points on the hyper-surface. The Liouville 
theory predicts the two point correlation function of the boundary primary operator Bfj[x) = 
ex.\){(54>{x)} to be 

1 



{Bp{Q)Bp{x)) 



|^|2A^ 



(18) 
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Fig. 10. Example of the 2-dimensional configuration with {N2, Ni) = (19, 17). 

where = /3(Q — Assuming the coordination number of a boundary vertex corresponds 
to be the boundary primary operator Bfj{x) with /? = 6 in the simphcial space, we measured 
the correlation function of coordination numbers of two vertices on the boundary using con- 
figurations simulated with n = = 0. We found that this correlation function showed a 
delta-function type singular behavior which had not been expected at all from the Liouville 
theory. The main reason of this property was thought to come from the singular nature of 
the surface created under the moves we have employed. The [AV, ASj-moves were designed to 
generate diagrams within connected and one-particle irreducible types, forbidding tad-poles 
and self-energy insertions at least for the manifolds with closed topologies. However, when we 
applied it to generate a surface with a boundary, it was found that the one-particle irreducible 
nature was lost. In fact, most of the configurations had many branches on the boundary such 
as the one shown in Fig. 10 . 

The matrix model also contains those diagrams which can be found in the table for the 
distribution of diagrams given in Appendix A. Under the condition of inflating space, most of 
the accepted moves are of the [AV, AS] = [1,1] type, and it obviously generates an uncorre- 
lated branches on the boundary. In order to keep diagrams within the one-particle irreducible 
type even for a surface with a boundary, we define a new set of moves as shown in Fig.ll, 
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Fig. 11. Typical [AV, AS']-moves for the Tutte algorithm. The vertices and a link on the boundary 
are denoted as w 




Geodesic Distance, D 



Fig. 12. The 2-point correlation function for the boundary operator B^"-\ 

consisting [AF, A5] = [±2,0], [1,1] and [—1,-1]. The configurations created by these moves 
coincide with the manifold studied by Tutte, which we introduce in Appendix B. The table 
for the distribution of diagrams generated by the Tutte algorithm clearly shows it prohibits 
branched surface. 

We measured the correlation function with the new configuration, but the results still 
exhibited the singular behavior. There seemed yet another important ingredient missing for 
obtaining the theoretically expected correlation function. We supposed it might be due to the 
improper definition of the primary operator in the simplicial space. In the previous simulations 
we chose the local coordination number as the simplicial primary operator. Instead, we defined 
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Fig. 13. 2-point correlation averaged 
over an ensemble of 10^ universes 
of typically N2 — 10"^ sizes. 
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Fig. 14. 2-point correlation function 
of super-long distances, choosing 
one universe at r = 5000, which 
has N2 - 10^. 



a new primary operator B^'^\i) in the simplicial space by the average of coordination numbers 
of those vertices covered vi^ithin n-steps from a vertex at i. Results are shown in Fig. 12 for n = 
1,2,3 using a configuration generated by the Tutte algorithm. The space dependence shows 
the predicted exponent, —2, beyond the separation 2n where the overlap of two operators 
terminates. We note that the correlation function using the old configuration remains singular 
even for the new operator. 

In quantum gravity the range of significant correlation is typically the Planck scale re- 
flecting the characteristic length of dynamical influence. When the space expands so rapidly, 
any dynamical perturbation cannot act on environments, and the quantum fluctuation will be 
expanded to macroscopic scale without much disturbance. In the construction of the partition 
function we have assumed that any possible configuration can be equally a member. As the 
consequence, the two point correlation function averaged over an ensemble of universes shows 
the short-distance nature of quantum mechanical scale shown in Fig. 13, where the long range 
correlation specific in each universe is wiped out after averaging over the ensemble consisting 
10^ universes of A''2 ~ 10^ sizes. 

In this respect one of the most interesting phenomena of the space-time geometry can 
be found in the observations of CMB anisotropies. An ultimate source of the anisotropies 
only we can imagine is the accidental creation of a seed universe by quantum fluctuation, 
and random expansion of space by accumulations of elementary units according to our model. 
When we measure the two-point correlation function, selecting one universe as COBE and 
WMAP observations did, it exhibits significant super-long distance correlation as shown in 
Fig. 14. The pattern of this fluctuation seems to be conserved in some degree during the 
expansion (Fig. 15) as we have imagined in analogy to a picture on an inflating balloon. We 
also calculate the angular power spectrum of the two point correlation function defined by 
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Fig. 15. Evolution of the 2-point function in the physical time t. 
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Fig. 16. Power spectrum at r = 4000, 5000 and 6000. 



|a;p with 

ai = j dOPiicos 9) f {cos 9), (19) 

where 9 = x/r{T) is the ratio of geodesic distance and the peripheral length at r. It shows 
the persistence of long distance correlation Fig. 16. 
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5. Discussions and Conclusion 

Stimulated by the reeent observations on the CMB anisotropies, a model to construct 
a quantized universe is proposed. Reexamining the DT-method developed in the simplicial 
quantum gravity, we extend the algorithm so that it can create a manifold with a boundary. 
Numerical simulation is made mostly on the 2-dimensional universe. Three types of universes, 
expanding, elongating and collapsing, are observed by varying two cosmological constants. 
Numerical results are closely examined by comparing to the matrix model and the Liouville 
theory predictions, which are considered to be the consistent theories of quantum gravity in 
2-dimension. The results exhibits fair agreements to the theoretical predictions. Physical time 
is defined based on the standard relation between volume and surface area, and the surface 
expands exponentially in the physical time. The exponential expansion is regarded as the 
quantum gravity realization of the inflation of the universe. Then, the two point correlation 
function on the equal-time hyper-surface is analyzed as the function of the geodesic distance. 
It shows the proper exponent expected by the Liouville field theory in the quantum scale 
after averaging over an ensemble of universes, while each universe has an extremely long dis- 
tance correlation. The anisotropies originated from quantum fluctuations and the subsequent 
inflationary expansion is one of the most significant findings of COBE and WMAP missions. 

Among many problems we need to clarify, the next step we should take is a numerical 
study on the universe in 4-dimension. We have already reported a few preliminary works on 
the algorithm and numerical results in several occasions. In 4-dimension there are eight 
[AV, AS'J-moves, corresponding to the four (p, g)-moves of the 3-dimensional DT and two 
types of attaching or removing a 4-simplex to or from the boundary for each (p, q)-move. 
Starting with a 4-simplex the Monte Carlo simulation is carried out for the partition function 
by constructing a Markov chain under the detailed balance condition accepting moves which 
fulfill the manifold conditions in exactly the same manner as the 2-dimensional simulation. 
The action to be used for the a priori probability pa — exp{—Sa) is assumed to be 

S = fiN4 + fi^Na + i^N2 + K^Ni, (20) 

where k and are referred to be the gravitational constant and the boundary gravitational 
constant. Extra two boundary terms are due to lack of relations existing in closed manifolds. 
Varying these parameters we can obtain several types of universes. For example, three kinds 
of universes similar to the 2-dimensional case are observed by varying /x and fi^ while two 
parameters, k and , are fixed to be zero. The expansion is exponential type in terms of the 
physical time in the simulation with fi = ji^ = 0. 

The equal-time hyper surface is considered to be a model for the very early universe. 
In order to appeal the fair possibility of the method we present here a comparison of the two 
point correlation function obtained from simulations and the WMAP data. The operator is 
chosen to be the scalar curvatures which are defined by the number of 4-simplices sharing one 
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Fig. 17. Comparison of 2-point correlation functions of the 4-dimensional simulation and the WMAP 
third year data. NormaUzation constant is /(O) for both cases. 



its boundary we create a B"^ cavity removing about half of the total 4-simplices covered 
within certain geodesic distance from an arbitrarily selected 4-simplex, and regard the 5^ 
boundary with the area as the last scattering surface(lss). Then, we choose one triangle i 
on the Iss and collect all the triangles {j{i,x)} with geodesic distance x from i measured in 
the dual lattice space. The two point correlation function is defined by 



where Ri and R are the scalar curvature at a triangle i and the average over N2 triangles, 
respectively. In order to compare to the observations, we also measure the two point correlation 
function of WMAP^ in the same manner as we have done for the simulation data: choose a 
point i, collect points x)} with the distance x apart from i, and calculate the temperature 
correlation function. The results are compared in Fig. 17. As far as the two point correlation 
function is concerned the agreement to the observation data is significant comparing to the 



Studies on the effects of extra two parameters to the structure of universe is one of the 
imminent problems. As in the 2-dimensional case, there is no reason to select specific values 
for these parameters. If the nature selects any, it should be the consequence of additional 
degrees of freedom. Therefore, the most important step we have to take is to go beyond the 
pure gravity by introducing the matter degrees of freedom. Without matter fields only the 
inflational expansion will continue forever. The universe is considered to have made a transition 
from the inflation to the big-bang by the deceleration effect of matter fields. Following the 



triangle. We use the simulation data of the space with the topology and A'4 ~ 10, 000. In 




(21) 



ACDM.2 
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same kind of procedures as the construction of simplicial manifold the additional variables 
are attached to each sub-simplex with appropriate rules for defining the neighbors. We should 
remember at this moment there exists no a priori laws of physics and thus no action function. 
The laws of physics which we expected to appear from the appropriate rules are the relativistic 
kinematics and quantum dynamics, which will be a primitive form of the relativistic quantum 
field theory. 

We have just started the study of the origin of a universe in 2-dimension as the first 
step. It is a small universe, but the necessary step toward the universe in 4-dimension. After 
all there should be universes in any dimension, if it is mathematically consistent. We recall 
the a mathematician's monologue^^ 'Is mathematics reality?'. We also agree his answer the 
reality lies outside us, and we physicists discover dynamical laws, which is a reflection of the 
mathematical reality.... 

We would like to give our sincere thanks to Dr.'s K. Hamada, A. Zamolodchikov, and 
H. Kawai for their stimulating discussions and especially N. D. Hari Dass for his careful 
reading of the manuscript.... 

Appendix A: Matrix model 

The one matrix model of cubic vertex is defined by the partition function. 




(A-l) 



with N X N Hermitian matrix M as the dynamical variable, where 



D[M] = Y[ dMii JJ d[^Mij]d[^Mij]. 



(A-2) 



i i>j 

Our main interest is one point Green's function, 




} 



(A-3) 



Fig. A-l. Correspondence diagram between the matrix model and the triangulated manifold 
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Table A-1. The coefficients /;,fe -the distribution of diagrams with Z-Unks and fc-triangles. 
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and its generating function F{j,g) = J^i^lid)^ ■ Expanding F{j,g) in powers of g and j 
^ k fiykd'^f coefHcient fi^k represents the number of diagrams with k vertices and I 
external Unes. The significance of this coefficient hes on the fact that the dual graph of a 
diagram corresponds to the two dimensional triangulated manifold with k triangles and I 
boundary links as Fig. A-1. 

After a few steps of mathematical manipulations for selecting connected planer diagrams 
excluding tadpoles and self-energy insertions, the equation for the generating function is ob- 
tained to be^^ 

gF^ -[g + j{l + s)]F + j{l + s) - gf +f = (A-4) 

with s{g) = gGs{g). The first few components of Green's function, Go{g) = l,Gi{g) = and, 
^2(5) = 1, reflect the constraints which we have imposed as the acceptable diagrams. This 
equation contains two unknown functions, F{j,g) and s{g). By requiring that they are both 
analytic at g = j = 0, s{g) is found to satisfy 

+ 3s^ + 8g^s'^ + 3s^ - 20g'^s + s+ 16/ - 5^ = 0. (A-5) 

Expanding these equations in terms of g and j recursively wc obtain the distribution fik. 
We list first several coefficients in Table A-1. From this table we notice that the branching 
polymer type diagrams such as the one with k vertices and k + 2 external lines are not 
forbidden. Diagrams we have accepted are the one particle irreducible type only for those 
without boundaries or a minimum boundary with 3-external lines. In the subsequent Appendix 
we introduce another type of generating function which counts only diagrams of one particle 
irreducible type even with a boundary. 

From the singularity analysis, s{g) is shown to be singular at g = gc = a/27/256, and 
F{j,g) becomes singular at g = gc, j = jc = \/3/4. Logarithm of these two critical couplings 
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correspond to the critical (lattice) cosmological constant Hc = 1.125 and the critical (lattice) 
boundary cosmological constant /i^ = 0.837. Expanding at the tri-critical point the first three 
leading singular terms of the generating function is given by 

F{j,9) ~ (j-Jc)-5(5-5c)^+2-^3-^(j-ic)-5(5-5c)^+2-'(i-jc)-^(5-5c)U.... (A-6) 

Using the binomial expansion formula, {x — Xc)°' = Yln^ni^c)^^ with ~ 
sin7ra ^_-[^-ji+a^-(i+a)p^2. + a){gc)" for a large n, the expansion coefficient of the distribution 
function can be written as 

1 1 

fn,m ~ m-2-V-^exp{-(// - //c)m} exp{-(//^ - //f )n}(l - 3 — + + •••)• (A-7) 

If we assume the exponential form as eq.(9) of the Liouville theory prediction the distribution 
function of the matrix model is expected to be 

fn,m ~ m-2-^n+°-^exp{-(Ai - /Xc)m} exp{-(^^ - //f )n}exp f-^— V (A-8) 
or in terms of the area A = ^c?m and the peripheral length I = an, it is written as 

Pm{1, A) ~ exp(-AA) exp(-A^O exp ( i-^^ . (A-9) 



Difference of the coefficient of the ^ from the Liouville theory by the factor 2 is also due to 
the symmetric factor. 

Appendix B: Tutte algorithm 

Tutte has been successful to obtain an equation for the generating function of the distri- 
bution of diagrams in the cubic vertex matrix model restricting within connected, and one- 
particle irreducible types. The generating function T{y,x) are expanded as n^rn,nx"y™' , 
where n is the number of loops and m-|-3 is the number of external lines. The exponent m and 
n are related to the number of vertices k and links I of the last Appendix ask = 2n + m + l 
and I = m + 3. The equation for T{y,x) satisfies. 



y'^T'^ + {x + xsy - y - y'^)T + y - xs = 0, (B-1) 

where the unknown function s{x) = T(0,x) is determined by the analyticity condition of 
T{y, x) dX X = y = 0, which leads 

(Sx^s^ + -20x5 + 32x - 1)2 + {8xs - i f = 0. (B-2) 

This equation is the same as Eq(74 -4), if we set s = xs and x = g^. The importance of 
the Tutte distribution function is the one-particle irreducible property of diagrams with a 
boundary. For example, the branching polymer type diagrams of m -|- 3 external lines with 
zero loops are forbidden as tjn,Q = 5m,o- The first few terms of the coefficients are listed in the 
Table B-1. 
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Table B-1. The coefficients tm,n -the distribution of diagrams with (m + 3)-Unks and n-loops. 
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The singularity of s{x) is identical to that of s{x) of Appendix A at a; = Xc = 27/256, 
and T(y,x) becomes singular at x = Xc, y = Vc = 3/16, or in terms of the variables g and j 
it is at Qc = jc = Logarithm of the surface critical couplings correspond to the 

critical (lattice) boundary cosmological constant /x^ = 0.5493, which differs from the value 
given in Appendix A. The generating function is expanded near the singularity similar to the 
last Appendix as 



11 , ,7 



T{y,x) ~ {j-j,)--2(g-g,)l+3-2i-f{j-jc)-H9-9c)-^+3\:^f{j-jc)--{g-gc)-^, (B-3) 
which lead the expansion coefficient of the distribution function to be 

/n,m ~ m-2-5n+°-^exp{-(// - //c)m} exp{-(/i^ - )n} exp(-^ — ). (B-4) 

The exponents coincide with the previous case, and the distribution function has the same 
asymptotic form as the previous case with the different value for fi^ and the coefficient in 
front of ^. 

We also can construct the moves which generate the diagrams satisfying the Tutte manifold 
conditions. They constitute four moves: [AF, AS'] = [2,0], and [1, 1], and their inverse moves, 
[—2,0], and [—1, —1]. They are shown graphically in Fig. 3. 

Appendix C: Boundary Liouville field theory 

The partition function of the Liouville field theory with the conformally invariant boundary 
is defined by^^ 

-^L = y -D[0] exp {-Abulk - Aound) ■ (C-1) 
The bulk Liouville action Aiyui^ is given by 

Abuik = ^J^ [g'^'^dadbct^ + QRcl) + 47rAe2^'^] ^/^(fx, (C-2) 
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where R is the scalar curvature associate with the background metric gab while Q{= b + 1/6) 
is the background charge. We consider only the geometry of a disk which can be represented 
as a simply connected domain T. The conformally invariant boundary action is written as 



■^bound _ 

IdT 



^[ [Qi^c^ + 27rA^eH5'/4de, (C-3) 
27r Jgr L J 



where K is the scalar curvature of the boundary associated with the background metric gab- 
The integration over ^ is taken along the boundary. 

Since the partition function is invariant under a constant shift + ^) the scaling 
relation, 

Z[X, \^] = e-^bXQz[Xe\ \^e''l\ (C-4) 

holds, where % is the Euler characteristic of the 2-dimensional surface. Choosing e*^ = A~^, 
the partition function is expected to be the homogeneous function, 

Z{\ \^\ = X^Q/^''Z[1, (C-5) 
V A 

This relation can be written in terms of the Laplace transformed variables defined by 

POO poo 

Z[X,X^]= dA dl e-^^-^''^Z[A,l], (C-6) 
Jo Jo 

where A and / are the variables corresponding to the area and the boundary length of the 
2-dimensional manifold, respectively. In terms of these variables the scaling relation reads, 

Z[A, I] = e-fbXQ-^zle-^A, e'^lH]. (C-7) 

Choosing e°' = A, we have 

Z[A,l]=A-i'^+i)z[l,^]. (C-S) 
According to Ref the scaling function is given by. 

For the case of pure gravity, b = a/2/3, with the disk topology, X = 1, the distribution function 
reads, 

Z[A,l] = A-'-'l-'-'exp (-^Q • (C-10) 
which differs from the matrix model prediction by a factor I. The partition function is given 

by 

Z[A, A^] = A^/^ ( 1 + V^h^^ J (l-^vW^^J, (C-11) 
apart from a numerical factor. 
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There are various correlation functions known analytically. Among them we write here 
necessary correlation functions for the primary vertex operators, Vb = e^^*^ and Bh = e^"^, which 
correspond to the bulk density operator and the boundary density operator, respectively, 

{B,{0)B,{x)) ^ (C.13) 

where the dimension = 1 for /3 = 6. The coefficients U{b) is related to the distribution 
function as U{b) ~ —§^Z[X,X^], 

/ \B \ 1/2 

U{b) - A^^ ( 1 + V^^^j , (C-14) 
while the coefficient d{b) is proportional to g^j^Z[A, A^], 

d{b) X^/^ + il + 2V^^^] , (C-15) 

within a numerical factor. 

Appendix D: Classical Liouville equation 

The Liouville field theory is considered to be valid only in the region where the stable 
vacuum exists. However, solutions of the classical equation of motion derived from the Liouville 
action, 

\dr,cl>f-{di:^f + A7rXe'''t>], (D-l) 

are known after J. Liouville. Here, the metric is taken to be Lorenzian and the sign of the 
cosmological constant is opposite to the Liouville field theory. The classical equation is written 
as, 

d^4> = dlcf) + ATrXbe'^'"!' (D-2) 
Changing the variables as y = ^ + r/ and z = ^ — r], the solution w{y, z) follows the equation 

dyd^w = -TrAte^''^. (D-3) 
The solution is known to have the following form, 

TTXb^{l- A{z)B{y)y' ^ ' 

where unknown functions A{z) and B(y) are fixed by the initial conditions. 
We consider a case when the initial conditions are given by 

e2'"^(o.«) = poiO, 5,</.(r7,e)|,=o = 0. (0-5) 

The second condition requires -B(^) = =byl(^). From the positivity of the density, e'^^'^^'^'^\ we 
need to choose the minus sign for > case. Writing p{ri, ^) = e^^'^^f^'^^ the solution is given 
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where 



cos2{i^(C-77)-F(C + ry)}' 



1/2, 



For the homogeneous case, po{Cj = Po (Const.), the solution is given by 



p{ri) 



Po 



u> 



(47rA62po)i/2. 
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(D.6) 
(D.7) 
(D.8) 



cos^(a;?7) ' 

Then the conformal scale factor a{ri) = e^'^ is written as 0(77) = cos'~^(a;77), and the physical 
time t defined by dt = a{r))dr] is expressed as 

1 + sin LOT] 



t = ±loglX^, (D.9) 
CO COS tori 

or cosh Lot = 1/ cos(a;?7). The boundary length a{r])lQ represented in terms of the physical time 
is expanding exponentially as l{t) = l{0) cosh(u;t). 
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